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We derive all single-field cosmologies with unit sound speed that generate scale invariant curvature 
perturbations on a dynamical attractor background. We identify three distinct phases: slow-roll 
inflation; a slowly contracting adiabatic ekpyrotic phase, described by a rapidly-varying equation of 
state; and a novel adiabatic ekpyrotic phase on a slowly expanding background. All of these yield 
identical power spectra. The degeneracy is broken at the 3-point level: unlike the nearly gaussian 
spectrum of slow-roll inflation, adiabatic ekpyrosis predicts large non-gaussianities on small scales. 



The observational evidence for primordial density 
perturbations with nearly scale invariant and gaussian 
statistics is compatible with the simplest inflationary 
scenarios. But is inflation unique? Are there dual cos- 
mologies with indistinguishable predictions? Such ques- 
tions are critical to our understanding of the very early 
universe. 

Inflation not only generates scale invariant and gaus- 
sian density perturbations, it does so on an attractor 
background. On super-horizon scales, the curvature 
perturbation on comoving hypcrsurfaces [1, 2], denoted 
by (, measures differences in the expansion history of 
distant Hubble patches [2]. In single-field inflation, ( 
approaches a constant at long wavelengths. In the strict 
k — > limit, ( — > 8a/ 'a, so the perturbation simply 
renormalizes the scale factor of the background solu- 
tion; such a perturbation can be removed by an appro- 
priate rescaling of global coordinates. For finite k, the 
perturbation cannot be completely removed, but differ- 
ent Hubble patches experience the same cosmological 
evolution, up to a shift of local time coordinates and 
a rescaling of local spatial coordinates. See [3] for a 
detailed discussion. 

Achieving both scale invariance and dynamical at- 
traction in alternative scenarios has proven challeng- 
ing. The £ equation of a contracting, matter-dominated 
universe is identical to that of inflation [4] , but ( grows 
outside the horizon, indicating an unstable background. 
The contracting phase in the original ekpyrotic sce- 
nario [5-10], with V{<j>) = -Vbe" 0/M , is an attrac- 
tor [11, 12], but the resulting spectrum is strongly 
blue [11-13]. A scale invariant spectrum can be ob- 
tained through entropy perturbations [14, 15], as in the 
New Ekpyrotic scenario [14], but this requires two scalar 
fields. 

The adiabatic ekpyrotic mechanism [16-20] proposed 
recently offers a counterexample: a single-field model 
for which the background is a dynamical attractor 
and generates a scale invariant (. The mechanism 



obtains for fairly simple potentials, such as V{4>) — 
Vb(l - e~*l M ) with V > and M < M Pl . Scale 
invariant perturbations are generated during the tran- 
sition when e = —H/H 2 = 3(1 + w)/2 rises rapidly 
from e <C 1, where the constant term dominates, to 
e w Mpj/2M 2 3> 1, where the negative exponential 
term dominates. 

Another counterexample proposed recently relies on 
a rapidly- varying, super luminal sound speed c s (t) [21- 
23]. See [24, 25] for earlier work. Even though the 
background is non-inflationary, £ is amplified because 
the sound horizon is shrinking. The growing mode is 
( — > constant, and the resulting 2-point function is scale 
invariant. 

The key lesson of these results is that relaxing some 
of the standard assumptions, such as w, c s w constant, 
opens up new possibilities for generating perturbations. 

In this paper, we derive the most general single-field 
cosmologies with c s = 1 that: i) yield a scale invariant 
power spectrum for (; and ii) are dynamical attractors, 
in the sense that £ — > constant is the growing mode solu- 
tion. These conditions imply a second-order differential 
equation for a(r) whose exact solutions we classify. 

The question of uniqueness is more than academic. If 
the Planck mission corroborates the predictions of the 
simplest single-field inflationary models, namely scale 
invariance and gaussianity, then the onus will be on 
theorists to determine whether inflation is unique in 
making these predictions. This work is an important 
first step in answering this critical and timely question. 

We find only three possibilities: inflation, with 
a(r) <~ 1/|t| and e « constant; the adiabatic ekpy- 
rotic phase [16, 18], with e ~ 1/r 2 on a slowly contract- 
ing background; and a novel adiabatic ekpyrotic phase 
on a background that first slowly expands, then slowly 
contracts [19]. At the 2-point level, therefore, the adi- 
abatic ekpyrotic phases are dual to inflation. The de- 
generacy is broken at the 3-point level: adiabatic ekpy- 
rosis generically predicts strongly scale dependent non- 



gaussianities, which limits the range of scale invariant 
modes that can be generated within the perturbative 
regime [18]. Thus, if Planck finds no deviations from 
gaussianity, our work will imply that any alternative 
theory must either invoke multiple degrees of freedom 
or use an altogether different mechanism to generate 
density perturbations. 

Any portion of these phases can be used to devise 
novel early-universe models. Such scenarios should 
explain the observed flatness and homogeneity, cither 
through inflation or through an ekpyrotic phase with 
e > 1 [5, 26]. Moreover, a reheating mechanism must 
be specified. In cases where the universe is contracting, 
the Null Energy Condition must be violated to bounce 
to an expanding phase, for instance within 4d effective 
theories [27]. 

For the purposes of this paper, however, we are solely 
interested in identifying all cosmological phases that 
generate, with a single degree of freedom, super-horizon 
perturbations compatible with observations — the can- 
didate duals to inflation. The idea of cosmological duals 
is not new [4, 11, 28], but we focus here on ( instead 
of the Newtonian potential [11, 28] and specialize to 
attractor solutions by demanding that ( — > constant. 



I. SET-UP 

Our starting point is the quadratic action for (, as- 
suming c s = 1: 
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pi 



drd 
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where z = ay/2e and primes denote derivatives with 
respect to conformal time r. This yields the mode 
function equation for the canonically-normalized vari- 
able v — z(: 



v'l + k 



Vk = , 



(2) 



where k is the comoving wavenumber. To generate a 
scale invariant spectrum from adiabatic initial condi- 
tions, it is sufficient for z to satisfy 
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r2 • 



Indeed, the solution to (2) in this case is 
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which implies that fc 3/2 |Cfc| = Vl + k 2 T 2 / ^2M ¥X z\t\. 
As t — »■ 0, fc 3 / 2 is independent of k, as desired. 

In addition to generating a scale invariant Cfcj our 
background must be a dynamical attractor. Since £fc <~ 
1/z|t| as k — > 0, the desired solution to (3) is 



V2 



(5) 



where m is an arbitrary scale. Combining (4) and (5) 
yields fc 3/2 |0c| = m\A + k 2 r 2 /2M Ph which is both 
scale invariant as r — > and constant as k — > 0. The 
observed amplitude of ( <~ 10~ 5 fixes m <~ 10 _5 Mpi. 

We pause to note that in an inflationary context the 
freeze-out or (-horizon |r| is usually identified with the 
comoving Hubble horizon, hr x = 1/aH — a/a 1 , but that 
more generally (e.g., when e varies rapidly) the Hubble 
horizon and the ^-horizon can differ greatly. 

Using the definition z = a\/2e, (5) implies 



(6) 



Moreover, we can rewrite e = —H/H 2 = dH 1 /At in 
terms of the comoving Hubble horizon h^ 1 = 1/aH as 

(h- 1 +r)' = e. (7) 

Combining (6) and (7) then gives a second-order differ- 
ential equation for a(r). Instead, we will cast these as 
a pair of coupled first-order equations. By differentiat- 
ing (6), 



(log^)' = -, 



-h. 



(8) 



Once we specify the signs of h and t, (7) and (8) become 
coupled ODEs for \h~ 1 \ and e. The behavior of (8) 
will depend strongly on the relative magnitude of the 
Hubble horizon \h~ 1 \ and the ^-horizon |r|. We will 
therefore say that the Hubble horizon is inside the (- 
horizon when < \t\, and outside the (-horizon 

when \h~ 1 \ > |t|. 

To solve these coupled equations, had an d esa must 
be specified at some fiducial time rgd < 0. To obtain a 
solution for a(r), we can set and = 1 by a spatial rescal- 
ing a — > Xa, r — > r/A. The equation of state is of course 
invariant, so e^d fixes Tad through (6). In practice, we 
will specify not |/ig d | but the ratio |/ig d |/|rfid|, which is 
invariant under the above spatial rescaling. 



II. SOLUTIONS 

While it is straightforward enough to integrate (7) 
and (8) numerically, as we have done, to guide our intu- 
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FIG. 1: Sketch of \h 1 \ for the contracting (dotted), expand- 
ing (dashed) and apex (thick dashed) branches of solutions. 



ition we also provide a series of simple, analytical argu- 
ments that explain the general features of the solutions. 
By varying over all possible initial conditions, we find 
three families of solutions, each of which is indexed by 
a single parameter and has finite duration, t\ < r < Tf . 
See Fig. 1 for a sketch of the solutions. 



A. Contracting Branch 

This case obtains if the universe is assumed contract- 
ing (/ifld < 0) at some fiducial time rgd < 0. Then, as 
long as h < and t < 0, (8) implies (logy/e)' > l/\r\, 
hence e increases monotonically. Meanwhile, (7) re- 
duces to = 1 — £) thus |ft _1 | increases when- 
ever e < 1, and decreases whenever e > 1. In fact, 
the bound from (8) implies that e must pass through 
e = 1, at which point \h ~ 1 \ hits a global maximum. 
A global maximum is a good point to specify a so- 
lution, so we denote the fiducial time in this case as 
T fid - > T max = —m^ 1 , where we set a max = 1 and 
£max = 1- All contracting solutions can therefore be 
indexed by the single parameter 



!^4=m|A-Ll>0. 



(9) 



Before 



7"max, a > 1, so e < l/m 2 T 2 ; after r max , a < 1, 
so e > l/m 2 r 2 . Integrating (7) therefore yields 



m\h~ l \ < f(r) < c, 



(10) 



where /(r) = c + 2 — m \t\ 
equalities saturated at r max 



- m |r| , with the in- 
Since m /i^ 1 < c, this 



implies that h cannot change sign for r < 0. Moreover, 
since f(r) vanishes at mr± = -(c+2± Vcv^+4)/2, 
h must diverge at finite r in both the past and the 
future of T max . Denoting the time of past and fu- 
ture divergences by r\ and T{, respectively, (10) implies 
t + < n < T max < Tf < r_ < 0. Over the interval 
n < r < Tf, a(r) contracts from oo to 0, so Tf marks 
a big crunch singularity; from (6), we conclude that e 
grows monotonically from to oo. 

The range of modes thus generated spans a factor of 
fcmax/fcmin = N/I r f| < l T +l/l r -l- From the definition 
of t± above, we have |r+|/|r_| = (mr + ) 2 < (c + 2) 2 , 
hence large values of c are required to generate a suffi- 
ciently broad range of scale invariant modes. From (9), 
this means that \h~ 1 \ must venture far outside the £- 
horizon, as sketched by the dotted line in Fig. 1. In this 
regime, e l/m 2 T 2 and a « 1, which is recognized as 
the adiabatic ekpyrotic phase proposed recently in [16]. 

Nearly all scale invariant modes are produced while 
\h~ 1 \ is outside the ^-horizon. Integrating (7) assuming 
e m l/m 2 T 2 gives ro|/i -1 | w /(r), or 



m\h 



For large c, horizon-equality (\h 
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(11) 



|r|) occurs at 



c 
2m ' 



1 
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(12) 



log |r< 



CC|-| 



hence this phase generates 
21ogc e-folds of modes. 

Because \h~ 1 \ is outside the (^-horizon during mode 
production, perturbations freeze out while inside the 
Hubble horizon and eventually exit Hubble by r cq _, 
when \h~ 1 \ re-enters the (-horizon. If a finite portion 
of this solution is used in a broader scenario, then some 
other dynamics must push these modes outside Hub- 
ble while maintaining scale invariance. In [16], this 
is achieved through an ekpyrotic scaling phase with 
e«c 2 /2>l. 



B. Expanding Branch 

Suppose instead that the universe is expanding 
(had > 0) at some fiducial time r^d < 0. It is helpful to 
rewrite (7) and (8) in terms of the gap A = \hr 1 1 — |r| 
between the Hubble horizon and the (-horizon. As long 
as ft. > and r < 0, (7) implies 



A' = e > . 



(13) 



Thus, when | h 1 1 is inside the (-horizon, corresponding 
to A < 0, the gap between the horizons narrows; when 
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\h~ 1 \ is outside the horizon, corresponding to A > 0, 
the gap between the horizons widens. Meanwhile, in 
this regime (8) becomes 

(logVe)' = \r\- 1 -(\r\+A)- 1 . (14) 

Unlike Case i), the evolution of e is no longer neces- 
sarily monotonic: when is inside the £- horizon, 
corresponding to A < 0, e decreases; when is out- 
side the C-horizon, corresponding to A > 0, e increases. 

It is straightforward to show that all solutions in this 
case must have emerged from a big bang singularity 
(where = 0) a finite time Tj < rgj in the past. 

In particular, is guaranteed to lie within the £- 

horizon at early times. Whether this remains the case 
subsequently depends on initial conditions. Qualita- 
tively, if \h~ 1 \ remains within the C-horizon, the solu- 
tion describes a universe that expands forever. This 
case, which includes the inflationary solution, is de- 
scribed below. If instead exits the C-horizon, the 
expansion inevitably comes to a halt at r = 0, and the 
universe enters a collapsing phase which terminates in 
a big crunch singularity. This apex solution is described 
in Case Hi). 

Let us now focus on the case where stays inside 
the C-horizon, i.e. A < 0. Since \h~ 1 \ < |r| < |rfid| 
for Tfid < t < 0, h cannot change sign as long as 
t < 0, hence a increases monotonically. From the dis- 
cussion below (14), e shrinks monotonically. In fact, 
since l^" 1 ! < |r| by assumption, \h~ 1 \ must hit zero at 
some Tf < 0. In other words, this case spans a finite 
time interval t\ < r < Tf, during which a(r) expands 
from to oo, while e shrinks from oo to 0. When e = 1, 
\h~ 1 \ reaches a global maximum, and, as in the con- 
tracting case, we can choose this as our fiducial time: 
Tfid -> T max = to -1 , where a max = 1 and e max = 1. The 
solutions can once again be indexed by c defined in (9) . 

Unlike the contracting case, c is bounded from above: 
|/imaxl nes inside the C-horizon by assumption, hence 
c < 1. For \h~ 1 \ to remain within the horizon subse- 
quently, we numerically find a tighter bound c < cq w 
0.52. As c approaches Co, Tf comes arbitrarily close to 
0. 

In fact, c w c is desirable to generate a broad range 
of modes, since fc max /fc m i n = | T"i | / 1 Tf j . In this limit, 
grazes the £- horizon, corresponding to e < 1 and |f?| = 
1 |dlne/di| <C 1. In other words, this case relies on 
a phase of slow-roll inflation to generate a broad range 
of modes. (Because we focus on exact scale invariance, 
this is a special case of slow-roll inflation. In particular, 
at linear order e and r] are related in such a way that 
n s — 1 = — 2e — 77 = 0.) The inflationary phase thus 



generates AT in f ~ log(l/m |rf|) e-folds of scale invariant 
modes, whereas mode production prior to the onset of 
the inflationary phase is negligible. Since Ih^ 1 ] < |r| 
throughout, modes exit Hubble before they freeze-out. 



C. Apex Branch 

In this case exits the C-horizon at some time 

Tcxit < after the universe emerges from the big bang 
singularity. Once this happens, there is no turning back 
— A becomes positive, and from (13) the gap keeps on 
growing for r < 0. 

From the discussion below (14), e attains a local 
minimum at horizon equality. The exit, defined by 
l^cxitl/l^xitl — 1) happens only once, so it is a conve- 
nient place to set a ox i t = 1. This family of solutions can 
therefore be indexed by a single parameter, e ox j t > 0. 

After horizon equality, the expansion inevitably 
comes to a halt at r = 0, at which time (the "apex") the 
universe enters a phase of contraction. The subsequent 
evolution can be deduced by noting that (7) and (8) 
are manifestly invariant under h — > —h, r — > — r. In 
other words, evolving forward in time when h > and 
t < is the same as evolving backwards in time when 
h < and r > 0. It follows that is guaranteed 

to reenter the ^-horizon, after which it will hit zero at 
finite Tf > 0, corresponding to a big crunch. 

To get a broad range of super-Hubble modes, we 
need e ox ;t <C 1 (corresponding to c w c ). This leads 
to a slow-roll inflationary phase, which occurs as be- 
fore while \h~ 1 \ grazes the ^-horizon, followed by an 
expanding adiabatic ekpyrotic phase [19], during which 
~> \t\, e <~ 1/t 2 and a(r) is slowly expanding. 
This solution thus includes two distinct phases of ap- 
preciable mode production. 

The inflationary phase generates iV; n f <~ e~ xit e-folds 
of scale invariant modes. Rescaling coordinates to set 
a = 1 when r = 0, outside the ^-horizon h~ x satisfies 

mhT x w -jttTV -1 - m(l - e oxit )T + e^ 2 . (15) 

Substituting in (7), we sec that the ekpyrotic phase 

with e w l/m 2 r 2 begins at T c k-bog ~ — m ~ le cx it 2 - This 
phase ends when Hubble re-enters the C-horizon, which 
from (15) occurs at r ok _ cnd w m _1 ^e cxit . 

The apex marks the end of mode generation. For 
t > 0, modes begin to re-enter the C-horizon, spoiling 
their scale invariance. Modes with fcr e k-cnd > 1 end 
up not scale invariant. The adiabatic ekpyrotic phase 
thus generates A" ok = log |r ck _bog|/|T c k-ond| ~ loge~ x - t 
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c-folds of scale invariant, super-Hubble modes. (Ar- 
bitrarily many c-folds can be obtained by ending this 
phase near r = while the modes remain within Hub- 
ble, but a subsequent phase would be necessary to push 
these modes outside Hubble while preserving their spec- 
trum [19].) 



III. NON-GAUSSIANITIES 

While the two non-inflationary branches which rely 
on a rapidly- varying e(i) yield power spectra identi- 
cal to that of inflation, the degeneracy with inflation 
is broken by non-gaussianities. The 3-point ampli- 
tude for the contracting adiabatic ekpyrotic mechanism 
was calculated in detail in [18]. The resulting non- 
gaussianities are strongly scale dependent and peak on 
small scales, with the dominant contribution growing 
as k . Since the 3-point calculation of [18] ignored the 
time-dependence of the scale factor, to a good approxi- 
mation the result applies equally well to the contracting 
or apex case. For completeness, we reproduce here the 
salient points of the 3-point amplitude calculation in 
the contracting case. 

To make contact with the results in [18], we intro- 
duce the parameter H = —m/c, where c was defined 
in (9). To see the physical significance of H , note 
that (11) implies that during the adiabatic ekpyrotic 
phase, —c/2m < t < — 2/mc, h^ 1 is within about a 
factor of two of its maximum value, /i^L 
follows that h~ x is nearly constant and 



H^ 1 . It 



h- 1 ~ H^ 1 (16) 

until near the very end of the phase. The parameter Hq 
is thus the characteristic Hubble parameter during this 
phase. Furthermore, the end points of the contracting 
adiabatic ekpyrotic phase, r cq+ to r oq _, are given by 



'eq+ 



(17) 



1 1 

2lf 5 T ° q - ~ c^Ho • 

Thus the long-wavelength cutoff for our calculations is 
T cq+ ~ Hq 1 . The short-wavelength cutoff is T e q_ , which 
is suppressed by a factor of 1/c 2 <C 1 relative to the 
long-wavelength scale. 

The cubic action for ( corresponding to a canonical 
scalar field with unit sound speed is given by, up to a 
field redefinition, [29] 

S 3 ~ ydM 3 x|e 2 CC 2 + e 2 C(VC) 2 -2eCVC-V X 

+ ^C 2 C + ^vc • v x v 2 x + ^v 2 c(v x ) 2 ], (is) 



where spatial derivatives are contracted with the Eu- 
clidean metric Sij, and X 1S defined as V 2 X = tQ. 
Moreover, following [18] we have ignored the time- 
dependence of the scale factor and set a ~ 1. At first 
order in perturbation theory and in the interaction pic- 
ture, the three-point function is 

( C(*,ki)C(t,k 2 )C(*,k 3 )> = 

-i r di / ([c(t,k 1 )c(t,k 2 )c(i,k 3 ), J ff int (t / )])(i9) 

J — oo 

where H mt = — L 3 , up to interactions that are higher- 
order in the number of fields, and to is chosen to be 
sufficiently late that all modes of interest have frozen 
out by this time. A natural choice in our case is 



1 



eq- 



c 2 H 



(20) 



As usual it is convenient to express the three-point func- 
tion by factoring out appropriate powers of the power 
spectrum and defining an amplitude A as follows 

(C(k 1 )C(k 2 )C(k 3 )) = (2^) 7 5 3 (5» > ( 21 ) 

where = fc 3 |£fc| 2 /27r 2 is the power spectrum for the 
curvature perturbation. 

The three-point function receives contributions from 
each interaction term in (18). The dominant contribu- 
tions, as shown in [18], are the last two terms in (18), 
both of which are 0(e 3 ). The next-to- leading contribu- 
tion is the f] term. We briefly review the calculation of 
these two contributions and refer the reader to [18] for 
further details. 

The e 3 terms give the combined interaction Hamilto- 
nian 



(22) 

Applying the canonical commutation relations, the 
three-point correlation function (19) in this case reduces 
to 

(C(k 1 )C(k 2 )C(k 3 )) £ 3 = *(2^) 3 <5 3 (5>)IIG^ ) 



dt 

kl dt 



dt 



+ perm. + c.c. > , 
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where the small imaginary part at t — > — oo projects 
onto the adiabatic vacuum state. Using the mode func- 
tions (4) and substituting e(t) ~ l/m 2 t 2 , it is easy to 
show that the integrand is a total derivative: 



d* 



3-iKt 



iKI 



AKt 



t 3 



- -r 6 H 3 

— c n , 



(23) 

where in the last step we have used (20) and taken the 
long wavelength limit K <C l^q-l^ 1 ss c 2 \H \, which is 
appropriate for the modes of interest. Putting every- 
thing together, the three-point amplitude is [18] 



K 2 
32H 2 



J^kf -J^kitf + 2fcifc 2 fc 3 ] • (24) 



As claimed, the 3-point amplitude is strongly scale de- 
pendent and peaks on small scales. 

The next-to-leading order contribution comes from 
the f] vertex in (18): 



"in, = / 'iV ' / / /s \'. 



(25) 



Using the fact that rj ~ —2m t 2 , the three-point am- 
plitude is given by, in the long wavelength (K <C c 2 \H \) 
limit, [18] 



j;*?-^^+tiij*3 ] • (26) 

i i^j 



This contribution scales as K/\H \ and is therefore sub- 
dominant relative to (24) on scales K > \H \. All 
other contributions to the three-point amplitude are 
suppressed by 1/c 2 <C 1 relative to (24). 

Following standard conventions, the three-point am- 
plitude translates into a value for f^ 1 ' , defined at the 




equilateral configuration: 



■I] 



equil. 
NL 



30 



A ki=K /3 5 K 2 



K 3 



144 H 2 ' 



(27) 



Unlike the power spectrum, the three-point function is 
thus strongly scale dependent: f^ 1 ' IS S C(i) on the 



largest scales (K <~ \Hq\) and grows as K 2 . Hence, as 
advocated, the degeneracy with inflation is badly bro- 
ken by non-Gaussianities. 

Since the perturbative parameter is /nlC> with ( <~ 
10~ 5 , perturbation theory breaks down for K > 
10 5 / 2 |ff |- In fact, on even smaller scales, K > 10 5 \H \, 
quantum corrections dominate the classical answer, sig- 
naling strong coupling [18]. As argued in [18], how- 
ever, these pathologies can be circumvented by modify- 
ing the e <~ 1/i 2 behavior before the dangerous modes 
are generated. In that case, the power spectrum for ( 
tilts strongly to the red and then flattens out at an ex- 
ponentially smaller amplitude with an acceptable non- 
gaussianity (/nlC 1) throughout. This results in a 
finite window (\H \ < K < 10 5 |7?o|) of scale invari- 
ant modes, which is sufficient to account for large scale 
structure and microwave background observations. 



IV. CONCLUDING REMARKS 

We have uncovered three distinct cosmological phases 
that yield a broad range of scale invariant modes: infla- 
tionary expansion, adiabatic ekpyrotic contraction [16], 
and adiabatic ekpyrotic expansion [19]. All three phases 
generate identical power spectra for £, and each is an 
attractor background. 

The degeneracy is broken at the 3-point level. 
The rapidly-varying equation of state characteristic of 
adiabatic ekpyrotic phases results in strongly scale- 
dependent non-gaussianities [18]. Our results imply 
that inflation is the unique single-field mechanism with 
unit sound speed capable of generating a broad range 
of scale invariant and gaussian modes. 

Forthcoming work [19] will extend the analysis to in- 
clude a general sound speed c s (t), the other degree of 
freedom at our disposal [25]. 
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